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5-h ' We find various lower and upper bounds for the index of Wente tori that contain 

a continuous family of planar principal curves. We then prove a result that gives an 
^S ' algorithm for computing the index sharply. 
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1 Introduction 

g : Hopf conned that eve ry constant m ean ™tn r e (CMC) dosed (olosed = compact 

J^, ■ to be either genus or embedded [H], but in general it is false, and the first counterex- 

amples (with genus 1) were found by Wente [We]. Abresch [A] then gave a more explicit 
representation for those Wente tori which contain a continuous family of planar principal 
curves, referred to here as symmetric Wente tori. Walter [Wa] later found an even more 
explicit representation for symmetric Wente tori, by noticing that if one family of principal 
curves are all planar, then the perpendicular principal curves each lie in a sphere. Finally, 
Spruck [Sp] showed that these symmetric Wente tori are exactly the same surfaces that 
Wente originally found. Pinkall and Sterling [PS] and Bobenko [B] went on to classify all 
closed CMC tori in M , and Kapouleas constructed closed CMC surfaces for every genus 
greater than one [Kl], [K2]. 

The index of minimal surfaces in R has been well studied. In this case there are no 
closed surfaces. Do Carmo and Peng [CP] showed that the only stable (index 0) complete 
minimal surface is a plane. Fischer-Colbrie [FC] showed that minimal surfaces have finite 
index if and only if they have finite total curvature. And for many minimal surfaces with 
finite total curvature the index has been found or has known bounds (for example, see [FC], 
[N], [MR], [T], [Cho]). 

However, less is known about the index of nonminimal CMC surfaces. It is known that 
the surface is stable if and only if it is a round sphere [BC] (analogous to the do Carmo- Peng 
result). It is also known that CMC surfaces without boundary in R have finite index if and 
only if they are compact [LR], [Si] (analogous to Fischer-Colbrie's result). However, other 
than the round sphere, there is no closed CMC surface with known index. (Unlike CMC 
surfaces in general, minimal surfaces have meromorphic Gauss maps when given conformal 
coordinates. Furthermore, for the index of nonminimal CMC surfaces we must consider a 
volume constraint that does not appear in the minimal case. This accounts for why more 
progress has been made in finding the index of minimal surfaces.) 

Our purpose here is to expand the collection of nonminimal CMC surfaces for which we 
have specific estimates of the index (to include surfaces other than just round spheres). We 
restrict ourselves to symmetric Wente tori, because they are the simplest known examples 



of nonspherical closed CMC surfaces in R , and can be represented nicely [A], [Wa], Hence 
they are the obvious candidates to consider in an initial investigation. 

In section 2, we briefly reiterate Walter's description of symmetric Wente tori, describe 
the Jacobi operator and the index problem, and discuss the eigenvalues and eigenfunctions 
of the Laplacian on a torus (which play an essential role in our discussion) . 

In section 3, we use Courant's nodal domain theorem to find initial lower bounds for the 
index, obtaining as a corollary that for every natural number N, there exist only finitely 
many symmetric Wente tori with index less than N. In section 4, we use upper and lower 
bounds of the Jacobi operator's potential function to determine upper and lower bounds 
for the index. 

In section 5, we use eigenfunctions of the Laplacian to create function spaces on which 
the Jacobi operator is negative definite, producing much sharper lower bounds for the index. 
These bounds imply that every symmetric Wente torus has index at least 7. Furthermore, 
the two simplest Wente tori (the first two surfaces shown in Figure 1) have index at least 9 
and 8, respectively. 

In section 6, we prove the correctness of an algorithm (a variation of the finite element 
method [FS]) for computing the index exactly, up to a possible change by 1 from the 
volume constraint. In this proof, we essentially show that on a flat torus, the spectrum 
of a Schrodinger operator A + V with C°° potential function V can be determined from 
the restrictions of the operator to function spaces spanned by a finite number of explicitly 
known eigenfunctions of the Laplacian A. 

Finally, in section 7, we implement the algorithm to give numerical estimates for the 
index of 17 different Wente tori (Table 3). These numerical estimates show that the lower 
bounds of section 5 are quite sharp, and show that it is reasonable to conjecture that every 
symmetric Wente torus has index at least 9. 

Acknowledgements. The author thanks S. Nayatani, H. Urakawa, K. Akutagawa, and 
M. Koiso for their helpful suggestions. 

2 Preliminaries 

2.1 Walter's representation [Wa] 

Consider a smooth conformal immersion X : J\A — ► R , where A4 is a compact 2-dimensional 
torus with the induced metric. If (x, y) are isothermal coordinates on M, then the first and 
second fundamental forms, Gaussian curvature, mean curvature, and Hopf differential are 

ds 2 = E(dx 2 + dy 2 ) , 11 = Ldx 2 + 2Mdxdy + Ndy 2 , 

The Gauss and Codazzi equations imply that <J? is holomorphic with respect to z = x + iy 
when H is constant. 

Since M is a torus, we may assume that M = C/T, where T is a lattice in the plane 
C. Since dz is a global 1-form on C/T, and since we will assume H is constant, we see that 
Qdz 2 is a constant multiple of dz 2 . 



We define a function F by 



HE = e F 



By a linear transformation of C, we may arrange that <&dz 2 = dz 2 , hence M = 0, L = e F +l, 
and TV = e — 1. Thus X has no umbilic points, and (x, y) become curvature line coordinates. 
Denoting the oriented unit normal by N, we may compute that 

11 ^11 



x w = -\f x Xx + \FyX y - (e F - 1)N , N x = H{1 + e~ F )X x , N y = H{1 - e~ F )X x 



2 2 



y ■ 



AF + 4i?sinhF = , (2.1) 

where A = -£-? + S-^. Thus finding the immersion X : C/T -» R is reduced to solving 
for F in equation (2.1). In the case of symmetric Wente tori, Walter determined F and the 
immersion X explicitly. In the process of doing so, he proved the following lemmas. 

Lemma 2.1 The set of all symmetric Wente tori are in a one-to-one correspondence with 
the set of reduced fractions ijn £ (1, 2). 

For each £/n, we call the corresponding symmetric Wente torus WV/ n . Each WV n has 
either one or two planar geodesic loops in the central symmetry plane: two loops if I is odd, 
and one loop if £ is even. Each loop can be partitioned into In congruent curve segments, 
and £ is the total winding order of the Gauss map along each loop. (The central planar 
geodesic loops are contained in the boundaries of the right-hand side figures in Figure 1.) 

Lemma 2.2 The function F for the surface W^/ n is 

F = 4arctanh(f \x)g{y)) , with 

f{x) = jcn k (ax) , g(y) = *fcnj,(ay) , 
where cn k (resp. cnj.) is the amplitudinus cosinus of Jacobi with modulus k (resp. k), and 

k = sin 9, k = sin 9 , for 9,9 € (0, ir/2) , and 9 + 9 < - , 



/ n l Ji L, T sin2fl / sin 20 

7 = Vtan#>0, 7 = vtan6>>0, a = « / AH — — , a = \IAH- 



sm2(9 + 9) \j sm2(9 + 9) 

There are two period problems for the immersion. One is a translational period prob- 
lem in the direction of the rotation axis of the surface, and the other is a rotational 
period problem around the rotation axis. Using elliptic function theory, Walter deter- 
mined that for all £/n, to solve the translational period problem of the surface, one needs 
9 = 65.354955354°. To then solve the rotational period problem, there is a unique choice 
of € (0°, 24.645044646°). The value of 9 (but not 9) will depend on £/n (see Table 2). 

Let X£ n (resp. y£ n ) be the length of the period of f(x) (resp. g(y))- 

Lemma 2.3 X : C/T — ► W^/ n is a conformal diffeomorphism , where 

r = span^{(nxg n ,0), (0, y£ n )} when £ is odd, and 

T = span%{(nx in /2, y in /2), (0, y in )} when £ is even. 
The curves { [xq ,y] \xq = constant} are mapped by X to planar curvature lines of We i n . 
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Figure 1: The surfaces W4/3, W3/2, and W 6 /5, with cut-aways to expose the inner parts. 
(Graphics by Katsunori Sato of Tokyo Institute of Technology.) 



The reason why T is different for £ odd or even is best explained in the introduction 
of [A]. (However, the notation "m/n" in [A] is only half the value of the notation "£/n" 
used here and in [Wa] . One can easily check that the characterization of T depends on the 
denominator's parity in [A]'s notation, and depends on the numerator's parity in [Wa]'s 
notation.) 

X maps the domain {(x, y) | < x < X£„/2,0 < y < ye n /2} to a fundamental piece of 
Wf/ n , bounded by four planar geodesies. The entire surface is composed of pieces congruent 
to the fundamental piece, in the sense that if one continues the fundamental piece by 
reflecting across planes containing planar geodesies, one arrives at the entire WV/ n . When 
£ is odd (resp. £ is even), W^/ n consists of the union of An (resp. 2ro) fundamental pieces. 

2.2 Index and the associated eigenvalue problem 

Let vt be the variation vector field at time t of a C°° variation M(t) C M of the surface 
Wi/ n so that M(0) = Wi/ n . Let A{t) be the induced area of M(t), and let 

u t := (v t ,N)-^3 , u:=u ■ 

It is natural to consider only those variations for which J(p, r UtdAt = id At is the induced 
area form on M(t), dA := dAo) for all t close to 0, corresponding to the fact that the volume 
inside a "soap bubble" is preserved by any physically natural variation of the bubble. Such 
a variation is volume preserving [BC]. So, the first variation formula implies that for volume 
preserving variations, 

A'(0) =H [ udA = 0. 

JC/T 

The second variation formula for volume preserving variations ([BC]) is 

A"(0) = f \Vu\ 2 + (2K - 4H 2 )u 2 dA = f u({-A + 2KE-4H 2 E)(u))dxdy, so 

JC/r JC/r 



A"(0) = / uC{u)dxdy , where (2.2) 

JC/r 

C(u) = -Au-V-u, V :=4HcoshF 

is the Jacobi operator with potential function V. (This operator is also computed in [PS], 
but with notation that differs by constant factors.) 

Note that in equation (2.2), we are integrating with respect to the flat metric on C/T, 
not with respect to the metric induced by the immersion. 



Definition 2.1 The index Ind(W(/ n ) of W^/ n is the maximum possible dimension of a 

i , v udA = and j<r\n 



subspace V C C°°(C/r) such that f(p, r udA = and J(p, r uCudxdy < for all nonzero 



ueV. 

By Lemma 2.4 of [BC], for every u satisfying J^.udA = 0, there exists a volume 

preserving variation of W^/ n with variation vector field uN on We/ n . Thus, loosely speaking, 
Ind(WV n ) is the maximum dimension of a space of variation vector fields for area-decreasing 
volume-preserving smooth variations, that is, a space on which A'(0) = and A"(0) < 
for all its nonzero variation vector fields. 



Let 1? = L 2 (C/r) denote the measurable functions u on C/T (or equivalently on W^/ n ) 
satisfying J(T, T u 2 dxdy < oo, and define the L 2 inner product (u, v) L 2 := /(P /r uvdxdy. If a 

sequence of functions Uj G L 2 converges strongly in the I? norm to a function u, we denote 
this by Ui ^ip. u as i — ► oo. Note that we define this L 2 space with respect to the flat 
metric on C/T, not the metric induced by the immersion. The following theorem is well 
known (see, for example, [U]): 

Theorem 2.1 The operator C = — A — V on C/r has a discrete spectrum of eigenvalues 

/3i < /3 2 < -.. / +00 
(each considered with multiplicity 1), and has corresponding eigenf unctions 

Vl ,v 2 ,... eC°°(C/r) 

which form an orthonormal basis for L 2 . Furthermore, we have the following variational 
characterization for the eigenvalues: 



Pj = inf sup / <pC(j)dxdy , (2-3) 

Vj V^ev,-,|H| i2 =i./C/r J 

where Vj runs through all j dimensional subspaces 0/C 00 (C/T). 

Lemma 2.4 If C has k negative eigenvalues, then Ind(W£/ n ) is either k or k — 1. 

Proof. Since Pk+i ^ 0, from characterization (2.3) with j = k + 1, we see there is no k + 1 
dimensional subspace of C°°(C/r) on which f(p . uCudxdy < for all nonzero functions u 
in the k + 1 dimensional subspace. Thus Ind(WW n )< k. 

Now let Vfc = spanj^i, ...,Vk}. Since (3j < for all j < k, Vk is a k dimensional space 
satisfying Jq. uCudxdy < for all nonzero u G Vfc. Define a linear map T : Vfc — > M. by 
F{u) = J(n /r udA. The kernel Kei(J r )= {u £ Vfc | J(n /r udA = 0} is a subspace of Vfc of 
dimension at least k — 1. Then, by definition, Ind(W^/ n )>dim(Ker(^ r ))> k — 1. □ 

Remark. It also follows from Lemma 2.4's proof that if there exists a j dimensional space 
Vj C C°°(C/T) with J£ /r uCudxdy < for all nonzero u £ Vj, then Ind(W^/ n )> j — 1. □ 

We define the nullity Null(£) of C to be the multiplicity of the eigenvalue of C 
Lemma 2.5 Null(C)> 6 for every We/ n - 

Proof. For any vector field tJonR generated by a rigid motion of R , the normal part 
(v, N)rn3N on YJii n is a Jacobi field on W^/ n , that is, C((v, N)^) = 0. Suppose Null(£)< 5. 

Then there is some such v entirely tangent to VVg/ n , since the rigid motions of R form a 6 
dimensional group. Choose such a v. For each p £ W^/ n , let 4> p {t) be the integral curve of 
v in W^/ n such that 4> p (0) = p. Since d(f) p (t)/dt = v^ u\, 4> p {t) is also an integral curve of 
v in R , hence W^ / n is invariant under the rigid motion that generates v. But W£/ ra is not 
invariant under any rigid motion, a contradiction, proving the lemma. □ 



2.3 Eigenvalues of the Laplacian 

By Lemma 2.4, our goal becomes to compute the number of negative eigenvalues of C. In 
order to do this, we use a convenient fact: we know explicitly the complete set of eigenvalues 
on and eigenfunctions m of —A. For C/T, with V = span^{(ai, 02), (61,62)}, the complete 
set of eigenvalues of —Aui = aiui are 

7 — 1 Tv> ((^262 -mia 2 ) 2 + (midi - m 2 bi) 2 ) , 

(ai02 — a26i) z v / 

with corresponding eigenfunctions 

c mi , m2 ■ (sin or cos) ( — 7-({m 2 b 2 - m 1 a 2 )x + (miOi - m 2 b 1 )y) ) , 

Vai6 2 -a 2 6i ) 

for mi,m2 G Z. We want these eigenfunctions to have L 2 -norm equal to 1, so we choose 
c mi , m2 = v/2/(|ai6 2 -a 2 6i|) if |mi| + |m 2 | > and c , = >/l/(|ai&2 -0261!). 

When £ is odd, by Lemma 2.3, we have a\ = nx£ n , a 2 = 61 = 0, 62 = ye n , and we 
can list these eigenvalues and eigenfunctions by choosing mi and m 2 in the following order: 
(ma.mi) = (0,0), (1,0), (0,1), (2,0), (1,1), (1,-1), (0,2), (3,0), (2,1), (2,-1), (1,2), 
(1,-2), (0,3), ... , always choosing sine first and cosine second. We can continue in this 
way to define all on and U{. The eigenfunctions U{ form an orthonormal basis for 1? . 

When £ is even, by Lemma 2.3, we have a\ = nxe n /2, a 2 = yin/2, 61 = 0, 62 = ye n , 
and again we can list the eigenvalues and eigenfunctions by choosing mi and m 2 so that 
(2rri2 — mi, mi) has the following order: (2m 2 — mi, mi) = (0,0), (2,0), (1,1), (1,-1), 
(0,2), (4,0), (3,1), (3,-1), (2,2), (2,-2), (1,3), (1,-3), (0,4), ... , always choosing sine 
first and cosine second. Again the eigenfunctions u\ will form an orthonormal basis for L 2 . 

Note that, with these orderings, we do not necessarily have on < a^ for i < i' . However, 
we still have lim^oo on = +00. 

Since we later refer to the ckj and Ui without restating their definitions, an unambiguous 
understanding of their ordering is needed, so for clarity we list the first 13 «j, U{ in Table 1. 

3 Application of Courant's nodal domain theorem 

In this section, we consider a Jacobi function on each surface W^/ n produced from rotation 
about the surface's central axis. Using the geometry of the surface, we can estimate the 
number of nodal domains of this function, and then the Courant nodal domain theorem 
gives the following lower bound for the index: 

Lemma 3.1 Ind(Wg/ n ) > 2n — 2 if £ is odd, and Ind(\Vi/ n ) > n — 2 if £ is even. 

Proof. Let v be the variational vector field of R associated to rotation about the axis of 
rotational symmetry of WV/ n . The normal part {v,N)-^3N is a Jacobi field on WV n , that 

is, v := (v,N)-rn3 satisfies C(v) = 0. Letting k be the number of negative eigenvalues of C, 
we may assume v = Vk+i- 

For odd £ (resp. even £), the set {[x, y] G C/T \ 2x/x£ n £ Z} is mapped by X to 2n (resp. 
n) closed planar geodesies in Wg/ n , each lying in a plane containing the axis of rotational 
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eigenfunctions 
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Table 1: The first 13 eigenvalues and eigenfunctions of — An« = ajUj. 

symmetry of WV n . Thus v^+i is zero along these In (resp. n) curves, so v^+i has at least 
In (resp. n) nodal domains. 

Courant's nodal domain theorem, which is valid in our setting (for example, the proof 
in [C] easily extends to a proof of this theorem for our operator C), says that the number 
of nodal domains of Vj is at most j for all j € Z . Choosing j = k + 1, we conclude that 
2n < k + 1 if £ is odd, and n < k + 1 if £ is even. Hence there are at least In — \ (resp. 
n — 1) negative eigenvalues of C if £ is odd (resp. even). Lemma 2.4 then implies the result. 
□ 

We immediately have the following corollary: 
Corollary 3.1 For every JVgZ, only finitely many W^/ n satisfy Ind(W^/ n )< N . 



We remark that Courant's nodal domain theorem can be similarly applied to produce 
lower bounds for the index of many of the surfaces described in [GBP], [HL], and [Kl]. 



4 Bounds for index derived from bounds for V 

As in the previous section, we again search for rough bounds for the index of the surfaces 

We/n- But instead of using a Jacobi function, this time we will use upper and lower bounds 

for the potential function, giving us both upper and lower bounds for the index. When the 

potential function is replaced by a constant upper or lower bound, we know the eigenvalues 

of the new resulting operator explicitly, leading to the lemma below. (This idea can be 

found in [FS], equation (15) of Chapter 6.) 

Let a p (i),a p (2)) ••• be the complete set of eigenvalues with multiplicity 1 of the operator 

— A on the flat torus C/T, defined as in section 2.3, but reordered by the permutation p{i) 

so that a p (i) < «p(2) < ••• /* +oo. (We use the reordering p(i) only in this section.) Define 

the constants 

Vmin = min V(x, y) , V max = max V(x, y) . 
[x,j,]e<C/r [x,j/]eC/r 

Consider the operators C m i n = — A — V m i n and C max = — A — V max . The complete set 
of eigenvalues with multiplicity 1 of C m in (resp. C max ) is then a p u\ = os p H) — V m in (resp. 

( ^p(i) a p(i) ~ ''max)- 

Recall that /% is a complete set of eigenvalues with multiplicity 1 of C. Comparing (2.3) 
with similar variational characterizations for the eigenvalues Q_ p u\ and a p u\ , we have 

°Lp{i) = Oi p {i) - V min > A > a p (i) = a p (i) - V max Vi . 
By Lemma 2.4, we conclude: 

Lemma 4.1 Choose \x € Z (resp. v G Z ) so that a p t^ < Vmin (resp. ol p i v ) < V max ) 
and otp^+i) > Vmin (resp. a p („ +1 ) > V max )- Then 

fi-l< Ind(W e/n ) < v . 

Lemma 4.1 gives explicit upper and lower bounds for the index, since we know the a p u\ 
explicitly (see Table 2). Since V max is generally much larger than V m i n , these bounds are 
very rough. However, when £/n is close to 1, the bounds get somewhat sharp, in the sense 
that the ratio of the lower and upper bounds is close to 1 (see, for example, the bounds for 
W21/20 an d W73/72 in Table 2). We produce much sharper estimates in the next sections. 

5 Specific spaces on which C is negative definite 

We now begin considering the Jacobi operator restricted to finite dimensional subspaces. 
The essential point is that if the operator restricted to some ./V-dimensional subspace is 
negative definite, then the index of Wi/ n is at least N — 1. We will use subspaces generated 
by a finite number of eigenfunctions of the Laplacian. This will lead us to stronger lower 
bounds for the index. 

Consider a finite subset {u± = u^ , ..., un = Ui N } of the eigenfunctions n« of — A on C/T 
defined in section 2.3 with corresponding eigenvalues dij = on ., j = 1, ...,N. Consider any 
u = YtiLidiUi G span{ui,...,ujv}, oi,...,ajv G K. 



/ uCudxdy = 2J o-i^j / Ui(—Auj — Vuj)dxdy = ^J ai(ctjSij — bij)a,j , 

J *i_-- /I ■ ■ -1 " ^s / L ■ ■ -1 



i,j=l I i,j=l 



where bij := J<r /T VuiUjdxdy. So Jjp ,„ uCudxdy < for all nonzero u £ span{ui, ..., ujv} if 

bij)ij=i ) ... ) N is negative definite. The remark following 



and only if the matrix Aat = {ctjdij 

Lemma 2.4 then implies the following result: 

Theorem 5.1 // the N x N matrix (aj5ij — &r/)ij=i,...,7v is negative definite, then 

Ind(W e/n ) >N-1. 

So, for example, consider the symmetric Wente torus W^- Choose N = 9 and choose 
u\ = ui, u 2 = u 2 , U3 = ""3, ua = u 4 , u 5 = u 5 , u 6 = u 7 , u 7 = u s , us = ug, and u 9 = uyj. 
Then, with H = 1/2, the matrix An is the following 9x9 negative definite matrix: 
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Hence, by Theorem 5.1, 

Ind(W 3/2 ) > 8 . 

This is the best we can do for the £/n = 3/2 case; any 10 x 10 matrix of the form (diSij — 
bij)i t j=i ! ,„ ! io will not be negative definite. 

For W4/3 we use N = 10 and the functions ui,...,uq,ui3 for ui,...,uio, producing the 
negative definite matrix (with H = 1/2) 



/ 



A 



N ■ 



V 



5.17 























-3.23 


\ 





-3.53 
































-3.53 
































-3.78 





-2.29 


























-3.78 





-2.29 




















-2.29 





-3.78 


























-2.29 





-3.78 
































-0.25 








3.23 























-2.21 
































-1.97 / 



Thus 



Ind(W 4/3 ) > 9 . 

For other surfaces we can use the same technique. Without showing the resulting ma- 
trices, we provide here some other examples: 

For W5/3 we can use N = 12 and the functions ui,U2,us,U5,...,ug,ui5,uiQ,uiT,U29- 



For W5/4 we can use N = 33 and ■ui,...,U23,'U27v>' u 35> u 45- 



For W7/4 we can use N = 16 and u±,U2 ,^3,^5,- ■■,ug,ui4,...,ui7,U27 ,^28 ^29)^45 ■ 
For W6/5 we can use iV = 20 and ui,...,u\g,U29- 
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Wi/n 





Xin 


Vln 


Vmin 


*max 


Lemma 
3.1's lower 
bound for 

Ind(HV) 


Lemma 
4.1's lower 
bound for 
Ind(W £/n ) 


Lemma 
4.1's upper 
bound for 
Ind(W £/n ) 


W3/2 


17.7324 


2.56 


4.21 


2 


123.447 


2 


2 


213 


W4/3 


12.7898 


3.28 


6.34 


2 


33.0184 


1 


6 


81 


W5/3 


21.4807 


1.76 


2.64 


2 


680.157 


4 


2 


743 


W5/4 


9.9285 


3.60 


7.90 


2 


17.9591 


6 


16 


169 


W7/4 


22.8449 


1.33 


1.94 


2 


2200.7 


6 


2 


1815 


W 6/5 


8.0983 


3.79 


9.18 


2 


12.4273 


3 


14 


83 


W7/5 


14.8978 


3.00 


5.38 


2 


54.5652 


8 


12 


351 


W 8/5 


20.1374 


2.08 


3.23 


2 


319.339 


3 


2 


433 


W 9 /5 


23.4867 


1.07 


1.54 


2 


5426.0 


8 


2 


3569 


W 7/6 


6.8332 


3.91 


10.3 


2 


9.65614 


10 


38 


189 


Wll/6 


23.8382 


0.89 


1.28 


2 


11308.6 


10 


2 


6191 


w 8/7 


5.9081 


3.99 


11.3 


2 


8.01729 


5 


24 


97 


W9/7 


11.1844 


3.47 


7.16 


2 


23.2959 


12 


28 


323 


Wio/7 


15.7491 


2.88 


5.02 


2 


68.2461 


5 


8 


277 


W11/7 


19.4966 


2.22 


3.50 


2 


238.928 


12 


6 


1037 


W12/7 


22.3044 


1.52 


2.24 


2 


1278.61 


5 


2 


1201 


W13/7 


24.0512 


0.77 


1.10 


2 


21012.8 


12 


2 


9863 


W21/20 


2.1359 


4.29 


20.2 


2 


3.54102 


38 


278 


491 


W73/72 


0.6005 


4.40 


39.1 


2 


2.3828 


142 


1962 


2353 



Table 2: Rough upper and lower bounds for Ind(W£/ n ). (x£ n , yi n , V m i n , and V max are 
computed using the value H = 1/2.) 

For W7/5 we can use N = 27 and ui,...,«ii,ui4,...,ui9,w 26 ,...,U3i,ti43, U44,ti45,u 6 5. 

For Wg/5 we can use N = 12 and u\,...,ut,u\q,.. .,^13,^29- 

For W9/5 we can use N = 20 and ui, 1^2,^3, ^5 ,...,ug,wi4v -,ui7, U26,- -,U29, ^43,^44,^45,^65- 

In each case, An will be negative definite. An is a diagonal matrix for W3/2, VV5/3, VV7/4, 
and W9/5, and is not diagonal for W4/3, W5/4, W6/5, VV7/5, and Ws/s- We conclude that 

Ind(W 5/3 ) > 11 , Ind(W 5/4 ) > 32 , Ind(W 7/4 ) > 15 , Ind(W 6/5 ) > 19 , 

Ind(W 7/5 ) > 26 , Ind(W 8/5 ) > 11 , Ind(W 9/5 ) > 19 . 

We could continue in this way for more WV n , but many computations are required to find 
sharp lower bounds. However, if we are interested in rougher lower bounds, then we may 
continue without so many computations. For example, for each of the surfaces Wg/7, W10/7, 
and W12/7, we can use iV = 9 and the functions u\, ...,U5,uio, ...,1*13 to create a negative 
definite diagonal matrix An, thus 

Ind(W 8/7 ) > 8 , Ind(Wi / 7 ) > 8 , Ind(Wi 2 / 7 ) > 8 . 
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From Lemma 3.1 and the above applications of Theorem 5.1, we conclude: 

Ind(W e/n )> 7 for all W e/n . 

Remark. In fact, numerical results in section 7 suggest that Ind(WV n )> 9 for all WV/ n . 
Since the symmetric Wente tori are possible candidates for minimizers of index amongst all 
closed unstable CMC surfaces, one may conjecture that, except for the round sphere, there 
does not exist a closed CMC surface in M. with index less than 9. □ 

6 A method for computing the index sharply 

In the previous section, we found finite dimensional subspaces on which the restricted 
operator is negative definite. Now, we will still use finite dimensional subspaces, but we will 
no longer confine ourselves to subspaces for which the operator becomes negative definite. 
Instead we will use the number of negative eigenvalues of the restricted operators on the 
subspaces to estimate the index of the surfaces from below. Again, as they are convenient to 
work with, we will use subspaces generated by eigenfunctions of the Laplacian. Furthermore, 
we will show that if we choose a sequence of subspaces whose union is dense in the full space, 
then the eigenvalues of the restricted operators converge to the eigenvalues of the original 
operator, providing us with a numerical algorithm for estimating the eigenvalues sharply. 

Let P m be orthogonal projection of L 2 to V m := span{«i, ...,u m }, the span of the 
first m eigenfunctions u, (as ordered in section 2.3) of —A. Hence P m {u = Y^£Li a i u i) 
= J2iLi a i u i = Ym=i( u -> u i) L 2U i- Note that P m Vj -^>i? Vj as m — > oo for any eigenfunction 



Vj of C 








Let C m = P m o C o P m . Then 








/ oo \ 


/ m \ m 






An 1 U = ^2 a i u i ) = P m° 


1 ^CliCUi J = ^2 a i(. a i$ij ~ 


Oij )Uj , 


where 


\ i=l / 


\i=l J i,j=l 






bij 


:= / VuiUjdxdy . 







So, on V m , C m is the linear transformation A m = (a,i5ij — hj)i<i,j,<m with respect to the 
basis {ui}^] =1 . The next lemma states that the linear map A m is the restriction of C to V m : 

Lemma 6.1 Let x = (xi,...,x m ),y = (yi,...,y m ) G K m . Let f = YT=i x i u i'& = YT=iVi u i- 
Then 

{x,A m y)-^rn = fCgdxdy and (x,y)^™ = fgdxdy . 



a 



Proof. j£ /r f£gdxdy = J27!j=i x iyj( a i^ij ~ hj) = {x,A m y)]^ m - Similarly, {x,y)j 
/C/r f&dxdy. 

Let A^ m < ... < Xm be the eigenvalues of A m with corresponding orthonormal eigen- 
vectors w\ m , ...,Wm . Note that since A m is a symmetric matrix, Xf 1 G R for all i. The 
following lemma is known ([RS], Theorem VIII. 3, or [FS], equation (23) in Chapter 6), but 
we include a brief proof here. 
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Lemma 6.2 A„ > A„ m > /3j for all j, m, m! such that j < m < m! . 
Proof. We have the following variational characterization: 

A; m) = sup inf / ipC m ipdxdy 

4>i,...,4>j-iev m \V>6V m n(span{</)i,...,^_i}) J -,||V|| i 2=wC/r J 

= sup inf / ip(P m oCoP m (tl;))dxdy) . 

^...^--isl 2 \V'GV m n(span{P m ( 9 ii),...,p m (^_ 1 )}) ± ,||v|| L 2=iJ(L/r / 

However, for -ip £ V m we have J(T/ r ^{P-m ° £■ ° Pm{^))dxdy = f(T, r tpCtpdxdy, so 

\j = sup inf / ipCipdxdy 

^...^igl 2 \^eV m n(span{^i,...,0 i _i}) ± ,||^|| L 2=i u 'C/r / 

So \j is nonincreasing as m increases. Also, (5j has the following characterization [U]: 

(5j = sup inf / ipCipdxdy 

</.i,...,^_iei 2 \V'ei 2 n(span{(}!>i,...,^_ 1 })i,||v|| L 2=iJC/r J 

Hence A y - > (3j for all m. □ 



We define 



7 7 - := lim A,- 



m^oo 



Am)-! 



This limit clearly exists, since {A^- }^" =7 - is nonincreasing and bounded below by (3j. 
Theorem 6.1 Jj = (3j for all j. 

Remark. The motivation for this theorem is that it gives us a method for estimating the 
eigenvalues (3j of C, since one can estimate jj numerically. Thus one can estimate the 
number of negative eigenvalues (3j < 0, and hence the index (see section 7). □ 

Before proving Theorem 6.1, we prove a crucial preliminary lemma, which essentially 
says that the Rayleigh quotient 

f£ /r (P m Vj)£(P m Vj)dxdy 

— '—, — -75- — of P m v j converges to the 

j£ /T {P m v j ydxdy 

j*C/r VjCvjdxdy 
Rayleigh quotient — -^- — -, — -75— — ; — = (3j of v 7 as m — ► 00, for all j . 

Ic/A v j) dxd y 

Lemma 6.3 f(n,(P m Vj)C(P m Vj)dxdy — ► f(n /r VjCvjdxdy = flj as m — ► 00 for all j. 
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Proof. (-Avj,Vj) L 2 = Ya?1=1 (-&((Vj,Ui) L 2U i ),(v j ,Ul) L 2Ul) L 2 = YX=i( v ji u i)h a i < °°- 
Similarly, {-AP m Vj, P m Vj) L 2 = J2"iLi( v j, u i)hai, hence 

„ m oo „ 

L (P m Vj)(-AP m Vj)dxdy = y2(v j ,Ui) 2 L 2ai^y2(v j ,Ui) 2 L 2a i = Vj(-Avj)dxdy 

JC/r i=1 . =1 JC/r 



as m — > oo. 



We complete the proof by showing /(P /r V(P m Vj) 2 dxdy — » J(P/ r Vv 2 dxdy as 



m — > oo. 



Let M+ = {[rr, j/] G C/r | (P mVj ) 2 -v 2 > 0}, and let M" = {[x,y] G C/r | (P m Vj) 2 -v] < 0}. 



C/r 



(P m Vj — Vj) 2 dxdy — ► 



{P m Vj — Vj) dxdy — ► and / (P m Vj — Vj) dxdy — ► 



as m — ► oo, so for M m equal to either M+ or M m , we have 



A/; 



\(P m Vj) 2 -v 2 \dxdy 



M; 



{{P m Vj) 2 -v 2 )dxdy 



m:, 



(P m Vj -Vj) 2 +2vj (P m Vj -Vj)dxdy 



< / {P m Vj - Vj) dxdy + 2\\vj\\ L2 ^ /r) \\P m Vj - ^|| L 2 ( C/r) ^ as m ^ oo 



u;, 



by the Cauchy-Schwarz inequality. Hence 



C/r 



V((P mVj ) 2 - vf)dxdy 



<V„. 



C/r 



\{P m Vj) 2 -v 2 \dxdy 



Vrn 



Mr, 



((P m Vj) 2 - v 2 )dxdy - V„. 



M~ 



((P m Vj) 2 - v 2 )dxdy 



as m — ► oo . 



□ 



Proof of Theorem 6.1. Since w\ m is a linear combination of Uj for j < m, we may con- 
sider w\ m as a function in C°°(C/r). Then J^.w^jCwudxdy = (w™,A m Wv — — 



S m ), 



(m) (m) 



(m) (m) y 



Xh Spv and /Q r ^ W cfedy = («;£ ; ,u^ ; )]R m = <W' b y 
so that -P m t>j = X^li a i7 "4 f° r an m anc ^ J' then 



Lemma 6.1. Choose a,-™ G 



C/r 



'P m Vj)C(P m Vj)dxdy = £(a^) 2 Aj m) - & 



i=i 



(6.4) 



by Lemma 6.3, and 



(PmVj, P m Vj) =J2(%] 



(™)\2 



1 



1=1 



(6.5) 



for all j G Z asm-> oo. 

Let -P™ be the orthogonal projection from L 2 to span{u^ , ...,w;> }. 

Choose ki,k 2 ,... such that ft = ... = (3 kl < (3 kl+1 = ... = (3 k2 < ft 2+1 = ... = ft 3 < 
ft 3 +i = ■■■ ■ Note that k\ = 1, since Courant's nodal domain theorem implies that any 
function in the eigenspace of the first eigenvalue has only one nodal domain, hence this 
eigenspace cannot contain two L 2 -orthogonal functions. 
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We will prove the result by induction on i, the subscript of ki, for i G Z . 
Base case: 71 = (3\, and P™vi —> L 2 v\ as m —> 00. 
We have 

m 

< (4^(71 - A) + (E(^V)(/3 2 " /3i) 

p=2 
mm m 

< (4?) 2 (7i - /3i) + ((E^?) 2 ) - m + i(E(«S } ) 2 ) - 11 • \Pi\ + (E(«i?) 2 )(/3 2 - Pi) 

p=l p=l p=2 

m m 

= (4?) 2 7i - a + (E(«!a) 2 )/?2 + i(E(«S ) ) 2 ) - 11 • i/3ii 

p=2 p=l 

m m 

< (E(«p?) 2 a^) - Pi + i(E(«S } ) 2 ) - 11 • ifti - 
P =i P =i 

asm^oo, by equations (6.4) and (6.5). Therefore, since /% — Pi > 0, we have 

m 

IIP" 1 -,, p „ II 2 — VV/.t'™)^ _. n 
W^i vi- P m vi\\ L 2 - 2_^{a p ,i ) ->■ U , 

p=2 

( a i 1 ) — * 1 ' an d (°i 1 ) (71 ~~ #l) — ^ as m — ► 00 . 

Thus 7! = /?!, and HPfH*! — v±\\ L 2 < ||P 1 m t'i — P m vi||^2 + ||P m t>i — v±\\ L 2 -> as m -> 00. 

Inductive step: If jj = fij and PJ^Vj —>l 2 v j as m ~> °° f° r an i < ki, then 
7? = ^7 an d P™ u i ~^L 2 v j as m ^ 00 for all j < fcj+i. 

Since P™fj —^2 Vj asfrnoo for all j < ki, we have P^fj — >£2 for all j £ [h + l, /Vt-i] 
asm->oo (since the Uj are orthonormal), so 

a H ~~ ¥ ^ ^ or an ' — ^* anc ^ an J ^ [^» + ^' ^*+i] • (6-6) 

Then for all j G [fej + 1, fej+i], we have 

fcj+l m 

°<( E («S ) ) 2 (7p-/?p)) + ( E («S ) ) 2 )(^ +1+ i-^ + i) 

p=fci+i p=fc i+ i+i 

< (E(«p?) 2 (7p - &)) + ( E (<#?) 2 )(/Wi " &h-i)+ 
p=i p=fc i+ i+i 

(E(«p?) 2 ^) + (( E («p™ ) ) 2 )-i)^ + i + (E(«S ) ) 2 i/3pi)+ 

p=l p=ki+l p=l 

m fcj+l 

i( E («S ) ) 2 )-i|-l^ + il = (E(«S ) ) 2 >)-^ + i+ 

p=fcj+i p=i 



( E («p7) 2 )^ +1+ i + (EK7) 2 I/ 3 pD+ 

p=fcl+l+l p=l 
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p=fci+i p=i 

(E^^l^D + K E («S ) ) 2 )- 1 I-I^+iH° 

p=l p=fci+l 

as m — ► oo, by (6.4), (6.5), and (6.6). Therefore, since (3k l+1 +i — flki+i > 0, we have 



E («^) 2 (7p-/3p)-0 (6.7) 



p=fci+i 



and 



|2 _ V^ (J m )\2 



IIP™ 7,- P 7,.|l 2 o - V^ Cn tmJN l 2 
p=k i+1 +l 

asm->oo. 

By (6.8) we have \\PjJ* ■ Vj - Vj\\ L 2 < \\Pj£ Vj - P m Vj\\ L 2 + \\P m Vj - Vj\\ L * — ► as 
m — ► oo, hence P™ i>j —^2 Uj asm-»oo for all j G [fcj + 1, fej+i]. Since, by the inductive 
assumption, PJ^Vj -^1,2 Vj asm^oo for all j < ki, it follows that also PjJ 1 Vj —^2 Vj as 
m — > 00 for all j < fcj. Hence P™ Vj —^2 Vj asm^oo for all j < /cj+i. 

Since P™ Uj ^^2 i>j as m — ► 00 for all j < hi + \, we have (P™ ^v^P™ v u )l2 — ► 

V as m ^ 00 for all M ,i/ < fe +1 . By (6.6), £^ +1 off off = (^V^i^,)^ " 

Sfii <4™ X™ -► £ M i/ as ra -► 00 when //, 1/ € [fcj + 1, fej+i]. So the matrix (a^ l J)k i +i<n,u<k i+1 
is almost orthogonal for large m. This implies, as (6.7) holds for all j £ [ki + 1, fcj+i], that 

7fc i+ i = ... = 7fc i+1 = /3fei+i- D 



7 Numerical computation of the index 

At first glance, the numerical computation of A m seems to involve a separate computation 
for each bij for 1 < i < j < m. Each bij is essentially a triple integral, since V is defined via 
the Jacobi en function, and en is defined via an integral. Numerically estimating m(m+ 1)/2 
triple integrals bij would require much computer time. However, most of the bij are 0, 
and there are clear relationships between many of the nonzero bij. So one wishes to find 
a minimal set of integrals that will determine all bij, hopefully having much less than 
m{m + l)/2 integrals, thus reducing the amount of required numerical computation. This 
is the purpose of this section. 
When £ is odd, we have 

I rnxtn r-yin 

6n = / / V(x,y)dydx , 

nX£ n yi n Jo JO 



and 



v^ t nxin [ yen rr, w. w 2vrax 2irby . , , 

b±j = / / V(x,y)(sm or cos)( 1 )dydx , 

nx£ n y en Jo Jo nx in y in 



for some a, b £ Z when j > 2, and 

2 f nx e« f ytn Tr , w ..2vrax 27r6y.. ..2vrcx 2vr(iy. , , 

Ojj = / / V{x,y){sm or cos)( 1 )(sm or cos)( 1 )dydx , 

nx in y£ n Jo Jo nx in y en nx in y in 
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Wi/n 


Theorem 5.1's 

lower bound 

for Ind(W £/n ) 


Theorem 6.1's 

estimate 
for Ind(>Ve /n ) 


4 


range of 

A m 's negative 

eigenvalues 


range of A m 's 

first 6 positive 

eigenvalues 


W 3 /2 


8 


10 or 11 


-4181 


(-35.4,-1.47) 


(0.059,0.65) 


W4/3 


9 


9 or 10 


A S i 


(-10.2,-0.63) 


(0.26,0.47) 


W5/3 


11 


11 or 12 


-4 8 5 


(-81.3,-50.6) 


(0.093,5.75) 


W5/4 


32 


33 or 34 


-4-145 


(-5.78,-0.035) 


(0.21,0.32) 


W 7/4 


15 


15 or 16 


-4-145 


(-193.7,-108.5) 


(3.88,11.3) 


W 6 /5 


19 


19 or 20 


-4 8 i 


(-4.37,-0.76) 


(0.12,0.19) 


W 7/5 


26 


26 or 27 


-4l45 


(-12.9,-0.72) 


(0.039,0.42) 


W 8 /5 


11 


11 or 12 


-4 8 i 


(-48.9,-0.89) 


(1.83,4.21) 


W 9/5 


19 


19 or 20 


-4l45 


(-294.3,-169.7) 


(9.38,21.1) 


W 7/6 





53 or 54 


-4i45 


(-3.82,-0.068) 


(0.11,0.37) 


Wn/6 





23 or 24 


-4l81 


(-445.3, -254.5) 


(10.4,21.9) 


w 8/7 


8 


34 or 35 


-4 8 i 


(-3.48, -0.024) 


(0.043,0.38) 


W9/7 


— 


53 or 54 


-4i45 


(-6.24, -0.074) 


(0.092,0.39) 


Wio/7 


8 


17 or 18 


-4 8 i 


(-12.1,-1.29) 


(0.33,0.83) 


Wn/7 


— 


34 or 35 


-4i45 


(-26.9, -0.037) 


(1.28,3.10) 


W12/7 


8 


13 or 14 


-4 8 i 


(-78.3,-57.4) 


(2.08, 10.6) 


W13/7 


— 


27 or 28 


-4l81 


(-503.0,-278.3) 


(24.8,37.3) 



Table 3: Sharper lower bounds and numerical estimates of Ind(W£/ n ). (The values in the 
final two columns are computed using H = 1/2.) 



for some a, b, c, d £ Z when i,j > 2. When £ is even, we have 

2 fnx ln /2 py ln 



>U 



and 



2^/2 



Kt 



nxenVen Jo 

nx ln /2 ryi n 



V(x,y)dydx , 



nXinHin Jo JO 

for some a, b £ Z with a + b even when j > 2, and 

rnx ln /2 ry e 



2irax 2-Kby 
V(x,y)(sm or cos)( 1 )dydx , 



nx ir> 



lllr, 



nxi n ye n Jo 



. 2irax 27rby . 2ircx 2-ndy 

V {x,y){sm or cos)( 1 J(sm or cos)( 1 )dydx , 

nxtn yen nx£ n ye n 



for some a, b, c, d £ Z with a+b and c+d even when i, j > 2. When £ is even, we have changed 
the domain of integration from {(x,y) \ < x < nxe n /2,y£ n x/nx£ n < y < y£ n x/nx£ n + ye n } 
to {(x, y) j < x < nx£ n /2, < y < ye n }, as this will not affect the value of the integral. 

We list the following facts without proof, as they follow easily from the symmetries 
V(x, y) = V(-x, y) = V(x, -y) = V((x e j2) -x,y) = V(x, (y e j2) - y) of V. 



If i + j is odd, then b 
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":n 



0. 



• If j is odd, and if a and b are not both even, then b 



1.7 



9 J1 



0. 
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• For bij / 0, 

kj = bji = Y^ c i,j,A,B • Z{A, B) , 

0<A,B<|o|+|6|+|c|+|d|,i4,B even 

where CijA,B G Z if i = j = 1 or i, j > 2, and Cij,A,B/V^ G Z if precisely one of i or 
j is 1, and 

1 f nx ^ fvtn 27rx A 2vry B 



I(A, S) := / / V(x, y)(cos ) A (cos —!L)"dxdy . 

nx in y in Jo Jo nx in y in 

The Ci t j t A,B are easily determined by repeated use of the identities cos (a ± (3) = 
cos(a) cos(/3) =Fsin(a) sin(/3), sin(a±/3) = cos(a) sin(/3)±cos(a) sin(/3), and sin 2 (a) = 
1 — cos 2 (a). 

So, for t odd (resp. £ even), if we compute A m=2 i2_2i + \ (resp. -4 m =4i 2 -4J+i) f° r some 
/ 6 Z , then we need compute only /(/ + l)/2 (resp. 2Z 2 — Z) integrals of the form T(A, B). 

Remark. We can compute the number of negative eigenvalues of A m . By Theorem 6.1 and 
Lemma 2.4 this number or one less than this number estimates Ind(WV/ n ). These estimates 
are listed in the third column of Table 3. The choice of A m is listed in this table's fourth 
column. For example, for VV3/2, we chose A m = A±gi, and computed numerically that 
Aisi has 11 negative eigenvalues. Thus Ind(W3/2) is estimated to be either 10 or 11, as 

written in the third column. And those 11 negative eigenvalues range from \\ ~ —35.4 
to A-q ~ —1.47 when H = 1/2, as written in the fifth column. (Changing H changes the 
eigenvalues by a factor of H, so the number of negative eigenvalues is independent of H.) 

The size of A m varies with the choice of £/n, since we choose A m as large as possible 
in each case without creating significant numerical errors. (Note that as the size of A m 
increases, the Cij t A,B quickly become very large, hence even very small numerical errors in 
the estimates of I(A, B) can result in large numerical errors, for large m.) 

By Lemma 2.5, Null(WV„)> 6, hence the first six nonnegative eigenvalues of A m will 
converge to zero as m — ► 00. This provides some indication of how close the eigenvalues 
of A m are to the first m eigenvalues of C Hence, in Table 3's final column, we include 
the range of the first six positive eigenvalues of A m . For example, the first six positive 
eigenvalues of the matrix Aisi for W3/2 range from A^ 2 ~ 0.059 to A^ 7 ~ 0.65 when 
H = 1/2. " ' □ 
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